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摘要 311P/PANSTARRS 是 一 颗 活 动 小 行星 , 具有 小 行星 和 花 星 的 双重 特征 , 是 中 国 “ 天 问 二 号 "的 探测 目 
标 之 一 . 311P/PANSTARRS 直 径 较 小 , 约 为 400 m, 非 引 力 效 应 可 能 会 对 其 长 期 动力 学 演化 产生 较 大 的 影响 . 
通过 假定 不 同 表面 组 分 , 研究 了 Yarkovsky 效 应 对 311P/PANSTARRS 轨 道 演化 的 影响 , 讨论 了 密 近 交汇 、 非 破 
坏 性 碰撞 和 YORP (Yarkovsky-O'Keefe-Radzievskii-Paddack) 效 应 等 非 引 力 效应 , 计算 了 小 行星 与 大 行星 密 近 
交汇 及 碰撞 概率 , 估计 了 311P/PANSTARRS 达 到 自转 周期 分 裂 极限 的 时 标 . 模拟 结果 显示 与 纯 引 力 模型 相 比 ， 
Yarkovsky 效 应 可 能 会 加 快 311P/PANSTARRS 离 开 当 前 共振 区 域 , 大 约 在 10 Myr 以 后 311P/PANSTARRS 会 
离开 当前 所 在 共振 带 , 在 表面 覆盖 风化 层 的 情况 下 有 机 会 通过 ve 长 期 共振 成 为 越 火 小 行星 ; 在 考虑 YORP 效 应 
的 情况 下 , 311P/PANSTARRS 在 2 Myr 时 标 内 可 达到 自转 周期 分 裂 极限 ; 在 考虑 Yarkovsky 效 应 及 YORP 效 应 
等 因素 的 情况 下 , 311P/PANSTARRS 在 10 Myr 时 标 内 仍 可 保持 其 动力 学 稳定 性 , 且 YORP 效 应 不 会 显著 影响 
其 半 长 径 偏 移 量 . 
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1 引言 RASEN, 因此 引入 了 “活动 小 行星 ”(Active As- 
茜 星 是 太阳 系 形成 时 期 洲 留 下 的 残 角 , 较 好 地 IDA ME AEST ee 
保留 了 太阳 系 形成 时 期 的 原始 信息 . SEMTE  CRABADE Ris DAE DARE 

OA EREN ee 质量 损失 的 机 制 有 很 多 种 , 比如 冰 的 升华 、 撞 击 喷 
在 轨道 动力 学 、 成 分 和 观测 特征 上 存在 一 定 的 过 发、 旋转 不 稳定 、 静 电力 、 热 分 裂 、 热 脱水 、 冲 让 
F, 传统 意义 上 认为 彗星 和 小 行星 是 两 种 不 同类 型 。 ”脱水 和 辐射 压力 等 四. 其 中 在 小 行星 主 带 因 冰 升 如 
的 天 体 . 但 近 些 年 的 观测 研究 发 现 , 在 小 行星 主 带 导致 周期 性 活动 的 小 天 体 称 为 主 带 彗星 (Main-belt 
中 , 许多 小 行星 也 会 出 现 类 似 彗 星 活 动 特 征 的 在 Comet), 而 因 其 他 因素 导致 短暂 活动 性 爆发 的 小 
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天 体 则 称 为 “破坏 小 行星 ”(Disrupted asteroid). X 
今 为 止 , 共 发 现 了 大 约 40 颗 活动 小 行星 . 其 中 大 部 
分 轨道 稳定 的 天 体 可 能 起 源 于 主 带 , 但 仍 存在 一 些 
天 体 轨 道 不 稳定 , 可 能 是 被 暂时 捕获 的 顽 星 . 因此 ， 


C 


土 0.09)，Hainaut 等 四 利用 3.56 m 新 技术 望远镜 
(New Technology Telescope, NTT) 2013 年 8 月 至 10 
月 的 观测 数据 以 及 S 型 小 行星 的 几何 反照 率 得 到 
311P 的 半径 上 限 范围 是 0.25-0.29 km. Jewitt 等 00 


对 活动 小 行星 轨道 演化 的 研究 可 探究 活动 小 行星 


在 2018 年 利用 哈 勃 望远镜 2015 年 3 月 至 7 月 的 数据 ， 


的 分 布 规律 , 从 而 为 搜寻 和 发 现 新 的 活动 小 行星 提 
供 理论 预测 ; 同时 , 通过 分 析 活 动 小 行星 的 轨道 迁 
移 机 制 , 为 推测 活动 性 背后 的 物理 原因 提供 限制 . 

311P/PANSTARRS, 又 称 P/2013 P5 (下 文 简 
称 311P), 2013 年 8 月 27 日 首次 由 Pan-STARRS (Pa- 
noramic Survey Telescope and Rapid Response 
System) 巡 天 计划 发 现 . 2013 年 9 月 ， 哈 勃 望远镜 
(Hubble Space Telescope，HST) 观 测 到 311P 呈 
DB DOR A ME BB AS, 这 一 活动 性 持 
续 了 约 5 个 月 (2013 年 4 月 15 日 至 2013 年 9 月 4 日 ). 根 
据 Jewitt 等 中 研究 表明 , 这 些 人 尘埃 尾 是 由 顽 核 快 
速 旋转 分 裂 产 生 的 . SERNA SAA IA, 
311P 的 Tisserand 参 数 T = 3.66, 在 小 行星 的 参数 
范围 内 (彗星 的 Tisserand 参 数 方 < 311). 这 些 特征 
都 表明 311P 是 一 颗 活动 小 行星 . 311P 是 “天 问 二 
号 ”的 探测 任务 之 一 ,“ 天 问 二 号 ”计划 于 2025 年 左 
右 搭 载 长 征 三 号 乙 运载 火箭 发 射 升 空 , 将 分 别 探测 
轨 近 地 小 行星 (469219) Kamo’oalewa 和 主 带 的 
活动 小 行星 311PLsl. 

311P 的 轨道 位 于 小 行星 主 带 的 内 缘 , 轨道 半 
长 径 约 为 2.189au, 有 着 中 等 偏心 率 ~0.142 和 较 小 
的 轨道 倾角 ~ 4.968°, 其 本 征 根 数 与 古老 的 Flora 族 
群 (~ 10° yrl 四 ) 小 行星 相似 , 可 能 是 其 中 的 成 员 之 
一 . Flora 族 群 是 一 个 具有 500 个 成 员 的 大 族群 , 通 
过 小 行星 光谱 检测 到 大 多 数 族 群 成 员 都 是 S 型 小 行 
EN, 另 一 方面 , 通过 Jewitt 等 回 对 311P 色 指数 的 测 
量 也 表明 311P 更 有 可 能 是 一 个 S 型 小 行星 . Jewitt 
等 g 利 用 哈 勃 望远镜 2013 年 9 月 10 日 和 23 日 的 观 
测 数据 计算 得 到 311P 的 绝对 星 等 (Hy = 18.69, 
18.54), 并 假设 其 几何 反照 率 与 Flora 族 群 的 平均 
反照 率 0.29 + 0.09 一 致 , 得 到 311P 的 有 效 半径 约 
为 (0.24 + 0.04) km. Jewitt 等 加 于 2015 年 利用 哈 勃 


试图 通过 光 变 曲线 数据 寻找 311P 快 速 旋 转 的 依 
据 , 但 在 此 数据 中 并 未 找到 相应 依据 , 只 得 到 自转 
周期 的 下 限 5.4 h. 同时 在 此 数据 中 得 到 的 绝对 星 
等 为 (19.14 土 0.02), 相应 的 有 效 半 径 约 为 (0.19 土 
0.03) km (假设 几何 反照 率 为 0.29). 

由 以 上 数据 可 以 看 出 , 相 比 于 其 他 活动 小 行 
星 ，311P 的 尺度 较 小 Lyapunov 时 标 较 短 (TI，= 
31.6 kyr), 这 意味 着 非 引 力 效 应 Yarkovsky 和 
YORP (Yarkovsky-O’Keefe-Radzievskii-Paddack) 
效应 对 其 轨道 的 影响 较 大 . 尽管 在 轨道 演化 中 通常 
认为 Yarkovsky 效 应 的 影响 较 小 , 但 近 些 年 通过 观 
测 发 现 小 尺度 小 行星 由 Yarkovsky 效 应 导致 半 长 负 
偏 移 是 非常 明显 的 Ha3. 因此 在 动力 学 演化 过 程 中 只 
有 同时 考虑 引力 和 非 引力 效应 , 才能 得 出 动力 学 稳 
定性 的 可 靠 结论 . 本 文 参考 以 上 数据 , 研究 311P 在 
Yarkovsky 效 应 影响 下 的 稳定 性 . 311P 的 许多 物理 
性 质 及 参数 尚 不 可 知 , 因此 可 以 利用 Fenucci 等 13 
最 近 使 用 的 方法 对 Yarkovsky 效 应 引起 的 半 长 轴 偏 
移 进行 统计 学 估计 . 内 容 安排 如 下 : 第 2 部 分 介绍 初 
值 设置 、 数 值 模拟 方法 和 理论 模型 依据 ; 第 3 部 分 
研究 稳定 性 、 可 能 的 轨迹 路 径 、 密 近 交 汇 以 及 在 
共振 带 的 停留 时 标 ; 第 4 部 分 讨论 了 密 近 交汇 、 非 
破坏 性 碰撞 以 及 YORP 效 应 对 轨道 的 影响 ; 第 5 部 
分 得 出 结论 . 


2 方法 
2.1 ” 纯 引 力 模型 和 初 值 设置 

本 文 所 采用 的 数值 积分 器 是 ChambersHg 开 
发 的 Mercury 程 序 包 中 的 BS (Bulirsch-Stoer) 和 
MVS (Mixed-Varieble Symplectic) 算 法 以 及 Fe- 
nucci 等 上 5 改进 的 Mercury 程 序 包 . 其 中 动力 学 模 


望远镜 2013 年 9 月 至 2014 年 2 月 的 高 分 辨 率 图 像 得 
到 311P 的 绝对 星 等 下 限 18.98 + 0.10, 相应 的 有 效 
半径 上 限 (0.20 + 0.02) km (假设 几何 反照 率 为 0.29 


型 的 引力 部 分 包括 八大 行星 以 及 月 球 ( 均 为 考虑 质 
量 的 大 天 体 , 月 球 考虑 为 地 月 系统 ), 311P 作 为 测试 
粒子 质量 为 0. 积分 器 中 设置 精度 为 10 17, 积分 步 


攻 
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长 设置 为 8d. 本 文 共 考 虑 了 两 种 动力 学 模型 一 纯 引 
力 模型 和 增加 Yarkovsky 效 应 的 引力 模型 (为 表述 
方便 下 文 有 时 简称 为 Yarkovsky 模 型 ). 在 纯 引 力 模 
型 下 , 测试 粒子 除 受 纯 引 力 影 响 外 , 其 他 非 引力 效 
应 均 被 忽略 . Yarkovsky 模 型 是 添加 了 Yarkovsky 效 
应 的 引力 模型 , 具体 详 见 2.2 节 . 


初 值 使 用 JPL (Jet Propulsion Laboratory) 的 
Horizons 所 给 根 数 , 所 有 大 天 体 取 颂 略 日 (Julian 


day) 2459000.5 JD (2020 年 5 月 31 日 ) 时 刻 的 瞬时 轨 
道 根 数 ,，311P 的 标 称 轨道 取 侍 略 日 2456894.5 JD 
(2014 年 8 月 25 日 ) 时 刻 的 瞬时 轨道 根 数 , 具体 数值 及 
误差 见 表 1, 表 中 o 为 轨道 半 长 径 、e 为 轨道 偏心 率 、 
i 为 轨道 倾角 、w 为 近日 点 幅 角 、Q 为 升 交 点 经 度 、 
M 为 平 近 点 角 . 为 了 对 观测 误差 范围 内 小 行星 可 能 
的 动力 学 特征 做 进一步 研究 , 将 lc 取 值 为 311P 的 观 
测 误差 . 轨道 测试 粒子 选取 方法 如 下 : 设 z* € ROW 
标 称 轨道 , T e 恨 6*6 为 对 应 的 协 方差 矩阵 , 该 矩阵 
来 源 于 JPL Small-Body Database. 则 标 称 轨道 的 
概率 密度 分 布 为 以 z* 为 均值 ， 以 万 为 方差 的 六 维 高 
斯 分 布 . 根据 高 斯 分 布 特性 , 若 w* c REEERE 
分 布 , 则 线性 形式 09 


Au* + 2%, (1) 


即 为 标 称 轨道 的 概率 密度 分 布 , 其 中 A © IRO* ry DL 
利用 Cholesky 分 解法 71 由 工 得 到 , HP = AAT, 其 
中 4 为 下 三 角 和 矩阵 . 根据 上 述 方法 , 首先 生成 标准 
正 态 分 布 随机 向 量 w*, 再 利用 (1) 式 可 得 到 六 维 高 
斯 分 布 的 轨道 测试 粒子 . 


2.2 ”Yarkovsky 效 应 及 初 值 设 置 

Yarkovsky 效 应 是 由 太阳 光 引 起 的 一 种 热 效 
应 , 可 以 影响 直径 小 于 30 km 的 小 行星 的 运动 轨 
迹 08-19. 旋转 天 体 吸 收 太 阳 热 辐射 , 然后 以 热 波 
式 重 新 辐射 出 去 , 从 而 产生 微小 但 持续 的 推力 . 1 
于 小 行星 的 热 惯 量 非 零 , 会 产生 辐射 延迟 , 因此 当 
天 体 旋转 时 推力 方向 与 日 心 位 置 矢量 不 平行 , 而 是 
略微 向 轨道 速度 矢量 方向 倾斜 , 进而 可 使 天 体 轨道 
半 长 径 发 生 偏 移 20-241. Yarkovsky 效 应 分 为 两 种 不 
同 效应 : 昼夜 效应 (da/dt)a 和 季节 效应 (da/d)s. 总 
的 偏 移 da/dt 为 两 种 效应 的 车 加 . 


INNS 


考虑 Yarkovsky 效 应 时 一 般 做 如 下 假设 : (i) 天 
体 为 球体 ，( 边 轨道 为 绕 日 圆 轨 道 ， (ii 天 体 表 面 满 
足 线性 边界 条 件 ,， (iv) 自 转轴 固定 . Yarkovsky 效 应 
的 大 小 取决 于 以 下 几 个 轨道 和 物理 参数 : a BE 
D、 密 度 p、 天 体 表面 导热 系数 KK 和 比热容 C、 自 
转轴 倾角 y、P、 吸 收 系 数 a 以 及 表面 发 射 率 e. 为 
了 估计 311P 的 半 长 径 偏 移 率 , 利用 公式 (具体 计算 
见 附录 ) 可 计算 311P 的 半 长 径 偏 移 率 [13: 
da da 
dt (SF 
其 中 大 部 分 的 参数 都 是 未 知 的 , 因此 我 们 根据 Fe- 
nucci 等 上 3 的 方法 , 利用 蒙特 卡 罗 模 拟 来 预测 (2) 式 
所 给 出 的 半 长 径 偏 移 的 可 能 值 . 


表 1 JPL Horizons 中 311P 在 2456894.5 JD 的 轨道 根 数 , 对 


应 的 1o 误 差 范 围 . 数据 来 源 : https://ssd.jpl.nasa.gov/. 
Table 1 The nominal osculating orbital elements 


) (@,D,p KO Pae), 2) 


of 311P and its corresponding lo uncertainties at 

epoch 2456894.5 JD, taken from JPL Horizons. 
Data sources: https://ssd.jpl.nasa.gov/. 
Orbital 


Rk Value lo uncertainty 

a/au 2.189085680217774 4.672 x 1078 

e 0.115558163253637 1.513 x 1077 
i/° 4.968041003605256 1.858 x 1075 
w/? 144.2928059076631 3.317 x 1074 
Q/° 279.2893239140268 2.781 x 1074 
M/°  39.87346296920205 1.677 x 1074 

2.2.1 Yarkovsky 效 应 的 初 值 设 置 


轨道 半 长 径 : Y 


长 径 的 1 不 确定 性 非常 小 , 大 


约 为 10-s au 的 数量 级 ( 见 表 1), 几乎 不 会 对 估计 半 


长 轴 漂 移 量 产生 相应 的 误差 , 因此 选用 标 称 轨道 的 


吻 切 根 数 . 
直径 和 密度 : 直径 可 以 由 以 下 公式 给 出 [9: 
b= ote, (3) 


其 中 pv 是 几何 反照 率 , 万 是 绝对 星 等 . 目前 研究 


未 能 给 出 准确 的 几何 反照 率 , 因此 本 文选 择 用 下 
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面 方 法 进行 估计 . 首先 对 311P 利 用 分 层 聚 类 分 析 


方法 (Hierarchical clustering method, HCM) 2324, 


振 之 间 , 与 木星 11:3 共 振 相 邻 , 因此 其 Lyapunov 时 
间 较 短 , ~ 31.6 kyr. Behrens 族 群 仅 有 4 个 较 大 小 行 


即 在 本 征 根 数 空间 中 确定 一 个 距离 闪 值 (分 层 


H3 > 
RK 


星 的 几何 反照 率 已 知 , 表 2 所 示 即 为 这 4 颗 小 行星 的 


方法 中 的 距离 , 单位 m/s), 选 出 在 距离 阔 值 范围 之 几何 反照 率 数据 . 以 此 4 个 小 行星 成 员 的 几何 反照 
内 的 所 有 小 行星 , 识别 其 所 在 的 小 行星 族群 . 本 征 率 为 准 计算 平均 值 , 得 到 Behrens 族 群 平均 几何 反 
根 数 是 一 定时 标 内 的 轨道 积分 得 到 的 平均 根 数 值 . 照 率 为 pv = 0.285 土 0.025, 并 设 为 311P 的 几何 反照 
Hsieh 等 叫 曾 利用 HCM 识 别 得 到 与 311P 相 关 的 所 率 . 与 Jewitt 等 B, 3, 1 所 假设 的 Flora 族 群 的 平均 反 
有 小 行星 , 并 命名 为 Behrens 族 群 . 本 文 利用 小 行星 照 率 相 比 , Behrens 族 群 可 能 更 为 年 轻 , 且 可 能 有 共 
数据 库 网 站 AstDys 再 次 对 311P 进 行 识 别 , 在 距离 同 起 源 , 因此 将 311P 的 几何 反照 率 设 为 Behrens 族 
阔 值 为 49 m/s 处 识别 得 到 25 个 成 员 的 Behrens 族 群 . 群 的 平均 反照 率 更 为 合理 . AI EH = 18.8 土 0.6 
图 1 是 Behrens 族 群 在 本 征 根 数 空间 的 分 布 情况 , I 


图 可 知 311P 位 于 6M-3S-10A 与 4MA1S-7A (“M 代 


至 于 311P 的 密度 , 由 


取 自 JPL. 再 利用 (3) 式 可 以 得 到 D = 433 + 121 m. 
于 对 311P 色 指数 测量 中 发 现 


表 火 星 Mars,“S” 代 表土 星 Saturn, 数字 表示 在 相同 其 更 偏向 于 S 型 小 行星 , 因此 密度 取 S 型 小 行星 的 平 
时 间 间 隔 内 完成 轨道 数 的 比例 ) 两 个 三 体 平 运动 共 均 密度 , 即 p = 2720 + 540 kg - m7325, 
= an 
! i 0.098 ¢ ia oe 
0.146} ! i A 。 noii 
1 i Bie esi 
l- i 1 2 lal 
0.145 + ! i 0.097 f TE 
1 1 1 | 
1 i dee me 
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| i Š I.e i 
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i | 0.095 + | 3 衣 
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0.142 F 1 i 4 1 党 党 
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FPE 6M-3S_10A 1 i 0.094 þe 6M-3s-10A 1 | 
pie AM41S-7A e l i | : | 二 4M+15-7A | | i : | 
2.178 2.182 2.186 2.190 2.194 2.178 2.182 2.186 2.190 2.194 
a, /au ay /au 
图 1 Behrens 族 群 在 本 征 根 数 空间 的 分 布 情况 , 左 图 是 本 征 半 长 轴 (ap) 和 本 征 偏 心率 (ep) 的 参数 空间 , 右 图 是 本 征 半 长 轴 和 本 征 倾角 正弦 


值 (sin 如) 的 参数 空间 
体 或 三 体 轨道 


园 


. 其 中 星 号 表示 311P 位 B., 
tHe NL, 虚线 表示 11J:3A, 点 线 表 示 6M-3S- 


点 表示 Behrens 族 群 成 员 , 点 的 大 小 与 成 员 的 有 效 


直径 大 小 


成 比例 . 标签 中 显示 的 是 与 大 行星 的 


OA, 点 虚线 表示 4M 十 1S-7A, 其 中 J 表 5 


RK EB, M 表 不 火星 ， S 表 示 土 性. 


Fig.1 The distribution in the proper element space of the Behrens family members. The left panel depicts the parameter space 


of semi-major axis (ap) and eccentricity (ep), while the right 


panel illustrates the parameter space of semi-major axis and 


sine of 


inclination angle (sin ip). Star represents the location of 311P, while circles represent the location of family members. Sizes of 


circle symbols are proportional to the estimated diameters of the corresponding objects. The labels indicate positions in 


resonance with the major planets. The dashed line represents 11J:3A, the dotted line represents 6M-35S-10A, and the 


dash-dotted line represents 4M+15S-7A, where J represents Jupiter, M represents Mars, and S represents Saturn. 


导热 系数 : 导热 系数 的 选取 依赖 于 311P 的 成 
分 、 孔 际 率 和 表面 温度 , 根据 不 同 的 表 塘 性 质 可 能 
会 变化 几 个 数量 级 , 不 确定 性 较 大 P9. 通常 情况 下 ， 
面 存 在 风化 层 时 , 导热 系数 较 低 , 一 般 在 0.0001- 


F 


0.1 W- m~. K-t 


围 内 ; 表面 为 裸露 岩石 情况 下 
导热 系数 较 大 , 一 般 在 1-10 W . m-1.K-! 范 围 
表面 含有 金属 成 分 的 情况 下 导热 系数 一 般 大 于 
10W-m!.K-!. 这 一 参数 通常 是 未 知 的 , 需要 


内 ; 


| 


额 
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外 的 数据 比如 红外 观测 或 Yarkovsky 效 应 探测 等 来 究 过 陨石 的 发 射 率 , 大 约 在 0.9-1 之 间 . 而 发 射 率 = 似 
估计 . 对 于 这 一 参数 的 选取 , 一 方面 考虑 到 311P 更 乎 并 不 会 显著 改变 半 长 径 偏 移 率 , 因此 选取 陨石 
倾向 于 是 岩石 成 分 的 S 型 小 行星 ; 另 一 方面 , Fenuc- ”发 射 率 的 平均 值 , 0.984. 吸收 系数 a 是 测量 小 行星 
ci 等 6 的 研究 表明 , 尺度 小 且 转 速 快 的 小 行星 可 能 。 ”表面 吸收 太阳 光 的 比例 , 由 a = 1 — 如 估计 式 给 
会 在 表面 保存 一 些 风 化 层 . 因此 本 文选 取 了 5 个 与 出 , A AAEREN, 一 般 而 言 小 行星 的 邦 
不 同 表面 成 分 相关 的 值 , 分 别 是 K © {0.001,0.01, ” 德 反照 率 小 于 0.15[29, 因此 a 大 约 在 0.85-1 之 间 . 同 
0.1,1,5} (W-m7!-K7"), 其 中 前 3 个 值 考虑 了 表面 样 , 由 于 其 他 参数 的 极 大 不 确定 性 , 在 本 模型 中 吸 
风化 层 存在 的 可 能 性 , 后 两 个 值 则 考虑 了 表面 为 不 收 系数 的 误差 不 会 对 半 长 径 偏 移 产生 较 大 影响 , 因 
同 孔 阶 率 的 裸 岩 时 的 可 和 ete. 此 将 吸收 系数 设 为 1. 

表 2 Behrens 族 群 已 知 几何 反照 率 的 小 行星 及 其 lc 误差 范 2.2.2 ”Yarkovsky 效 应 估计 
o B |e eae 利用 上 一 小 节 所 描述 的 Yarkovsky 模 型 , 生成 
able e known geometric albedo of the 、 、 
Behrens family and fac corresponding lo errors ed i Mi Ap EY 3 数据 ， UF ls gee ia 
and sizes. Data sources: 311P 的 六 长 偏 移 率 da/dt. 然后 使 用 核 密度 估计 
ey age 得 到 不 同 导热 系数 到 值 的 连续 概率 密度 函数 , 结果 

Asteroid lo Dia 如 图 2 所 未， 7 p 

name uncertainties 对 于 前 4 个 小 于 等 于 1 的 导热 系数 有 两 个 可 能 
1651 0.381 0.052 8963 性 较 大 的 值 ， 个 正 值 和 一 个 负 值 . 这 是 | 于 (2) 式 
Se a sag 中 昼夜 效应 项 占 主 导 ， 4K = 
0.001 W- m7! - KHF, 正 负 两 侧 关 于 da/dt = 
13410 0.237 0.043 2654 0a /Mys EAE 的 ， 随 着 KK 值 的 增 大 ， 负 值 占 
88630 0.316 0.161 1240 ‘ee 


比热容 : 比热容 同样 取决 于 小 行星 的 成 分 和 表 
面 温度 , 但 其 变化 范围 并 没有 导热 系数 广 . 对 于 表 
面 履 盖 风 化 层 或 裸露 的 岩石 主 带 小 行星 ， 该 值 一 
般 取 600-700 J- kg! . K-27), 另外 , 与 导热 系数 相 
比 , 比热容 的 不 确定 性 对 于 da/dt 的 影响 较 小 . 

自转 轴 倾 角 : 311P 的 自转 轴 倾 角 是 未 知 的 , 因 
此 选取 cos y € [一 1, 1] 的 随机 值 . 

自转 周期 : 自转 周期 一 般 
然而 由 于 观测 时 间 有 限 , 目前 
光 变 曲线 数据 直接 证 明 311P 为 快速 旋转 的 小 行星 . 
因此 对 于 自转 周期 的 估计 , 假设 赤道 速度 为 逃逸 
速度 , 利用 ve。 = V2Gm/r, 其 中 G 是 万 有 引力 常 
数 , r 是 311P 的 半径 , m 是 311P 通 过 密度 (p = 2720 
+ 540 kg.m-3) 计 算 的 估计 质量 , 可 得 到 自转 周期 
的 临界 值 Powit ~ 1.4+0.2 h. 因为 311P 是 一 颗 旋 转 
裂变 的 活动 小 行星 , 可 将 自转 周期 初 值 取 其 临界 
值 . 

表面 发 射 率 和 了 吸收 系数 : Ostrowski 等 23] 4h 


由 光 变 曲线 来 确定 ， 
研究 中 没有 相应 的 


三 | 
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逐渐 增 大 , MACK = 5.0 W:m-1.K-1 时 几乎 全 
侧 仅 集中 在 非常 小 的 da/dt (~ 
这 种 情况 的 出 现 主要 是 
应 项 相当 , 而 季节 


小 且 负 值 


侧 


出 现 两 个 峰值 , 正 值 


0.0001 au/Myr)4b 
季节 效应 项 与 入 
导致 半 长 径 


HTA 
效应 项 总 


最 终 造成 半 


Kil 


RE 
局 移 量 偏向 负 值 的 现象 . 对 于 小 


热 系数 ( 即 K = 0.001, 0.01 W- m~t. K- 
半 长 轴 偏 移 较 快 ， 人 
FK =0.1W-m"!-K—!, KAW 
对 于 较 大 导热 系数 KK = 
要 集中 在 负 值 侧 ， 约 为 -0.0002 au/Myr 以 
K-:! 时 两 个 负 峰 值 ， 约 


则 主 


RK = 50W .m-! 


1.0W .m-! 


为 一 0.0004 au/Myr 和 一 0.0001 au/Myr. 


3 结果 


利用 2.1 


子 加 入 2.2 节 吕 
进行 轨道 演 


PP 所 描述 的 Yarkovsky 效 应 初始 参数 


化 模拟 . 


夜 效 
向 内 偏 移 ， 


Ey 


导 
1) 的 情况 ， 
0.0008 au/Myr， 而 对 
上 0.0003 au/Myr. 
K-t, 


节 中 的 方法 生成 500 个 测试 粒子 在 纯 引 


力 模型 下 进行 轨道 演化 模拟 , 然后 用 同样 的 测试 粒 


时 
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Probability Density Function 
Probability Density Function 
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图 2 ”模拟 增加 不 同 导热 系数 的 Yarkovsky 效 应 后 引起 的 311P 半 长 轴 偏 移 率 da/dt, 其 中 红色 曲线 表示 利用 核 密度 估计 拟 合 的 连续 概率 密度 函数 . 


Fig.2 Modeled semi-major axis drift rates da/dt caused by the Yarkovsky effect for asteroid 311P, provided for various values 


of the thermal conductivity K. Red curves represent continuous probability density function using the kernel density estimation. 
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ti 
Eb 
$ 
m 
= 
$ 


行星 311P/PANSTARRS 的 动力 学 特性 分 析 2A 


3.1 ”长 期 演化 

纯 引 力 模 型 下 的 标 称 轨道 如 图 3 所 示 . 6M- 
3S-10A 的 共振 中 心 位 于 a ~ 2.1886 au, 共振 宽度 
约 为 0.0012 au，4M+1S-7A 的 共振 中 心 位 于 a ~ 
2.1895 au, 共振 宽度 约 为 0.0017 au, 因此 311P 位 
于 [2.1874，2.1912] au 的 共振 区 域 中 . 从 图 中 可 知 
在 ~ 17 Myr 时 ，311P 从 6M-3S-10A 三 体 共振 跳跃 
至 4M+1S-7A 三 体 共 振 , 之 后 在 两 个 弱 三 体 共 振 间 
来 回 跳跃 , 但 311P 在 两 个 弱 三 体 共振 的 影响 下 混沌 
扩散 并 不 明显 , 偏心 率 、 倾 角 与 近日 距 g 基 本 保持 
稳定 , 分 别 在 e € [0.05, 0.22]、 i € [2°, 8°]. q € [1.70, 
2.05] au 范围 内 , 这 说 明 两 个 弱 三 体 共振 并 没有 
重合 而 导致 混沌 扩散 . 图 4 是 纯 引 力 模 型 中 500 
个 测试 粒子 在 (-1, 1) Myr 时 间 区 间 内 轨道 半 长 


径 的 演化 范围 , 图 中 黑色 曲线 是 所 有 粒 子 半 长 径 
的 平均 值 , 灰色 区 域 是 所 有 粒子 半 长 径 的 标准 差 
范围 , 偏心 率 和 倾角 在 (一 1, 1) PO 
明显 . 起 始 时 间 附 近 偏 差 极 小 , 说 明 所 有 测试 粒 
子 都 有 相似 的 轨道 . 在 经 过 大 约 200 kyr 的 演化 后 ， 
标准 差 开 始 逐 渐 增 大 , 但 平均 值 仍 在 2.1886 au 附 
VE, 这 说 明 在 纯 引 力 模 型 的 情况 下 受 6M-3S-10A 共 
振 的 影响 更 大 . 图 5 则 是 增加 Yarkovsky 效 应 后 测试 
粒子 的 轨道 演化 范围 . 通过 对 比 , 在 增加 Yarkovsky 
效应 后 测试 粒子 轨道 在 大 约 60 kyr 时 标准 差 就 开始 
BOK, 平均 值 略 高 于 纯 引 力 模 型 下 测试 粒子 的 轨 
道 . Yarkovsky 效 应 的 引入 导致 半 长 轴 的 平均 值 有 
~ 2 x 1073 au 的 偏差 , 而 偏心 率 与 倾角 在 此 时 间 尺 
度 内 几乎 没有 变化 . 
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93 纯 引 力 模型 下 ,311P 标 称 轨道 在 50 Myr 的 演化 . 从 上 至 下 分 别 为 半 长 径 、 偏 心率 、 倾 角 以 及 近日 距 . 


Fig.3 The long-term evolution of 311P’s nominal orbit over 50 Myr, computed within a purely gravitational model. From top 


to bottom, the panels depict the semi-major axis, eccentricity, inclination, and perihelion distance. 


为 了 研究 Yarkovsky 效 应 对 311P 长 期 演化 的 
影响 , 我 们 分 别 绘制 了 纯 引 力 模 型 与 考虑 K = 
0.001 W- m-1.K-! 的 Yarkovsky 效 应 后 100 个 测试 
粒子 的 轨道 分 布 R(a,e, 让 . 首先 将 轨道 根 数 空间 
(ae 让 分 为 Aa x Ae x Ai, 其 中 Aa = 0.004 au, 
Ae = 0.001, Ai = 0.1°, 然后 对 测试 粒子 进行 轨道 


演化 , 演化 时 标 为 50 Myr, 输出 时 间 间 隔 为 10 kyr, 
并 计算 落 在 每 个 区 间 内 的 次 数 , 最 后 标准 化 该 分 布 
使 得 


i R(a,e,i)dadedi = 1. (4) 
实际 上 该 轨道 分 布 给 出 了 测试 粒子 在 每 个 给 定 的 
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轨道 根 数 空间 所 停留 的 时 间 间 隔 . 在 纯 引 力 模型 下 
以 及 Yarkovsky 模 型 下 的 轨道 分 布 如 图 6 和 图 7. 两 
种 模型 最 密集 的 区 域 都 位 于 a ~ 2.1886 au 附 近 , 偏 
心率 小 于 0.2、 倾 角 小 于 8°. 当 Yarkovsky 效 应 引入 
后 , 该 区 域 的 分 布 密度 比 纯 引 力 模 型 下 略 小 . 在 这 

一 区 域 以 外 , 测试 粒子 有 10% 分 散 到 近 地 轨 道 区 域 ， 
甚至 与 金星 密 近 交汇 . 


—1000 -800 -600 -400 -200 0 200 400 600 800 1000 


t/kyr 
图 4 纯 引 力 模型 下 , 311P 的 500 个 测试 粒子 在 (一 1, 1) Myr 内 的 演 
化 . 其 中 黑色 曲线 代表 所 有 测试 粒子 半 长 轴 的 平均 值 , 灰色 区 域 则 是 对 
应 半 长 轴 的 标准 差 . 
Fig.4 Evolution of 500 orbital clones of 311P in the time 


interval (—1, 1) Myr, based on a purely gravitational model. 


The central black line represents the average of the 
semi-major axis and the gray filled area represents the 


corresponding standard deviation. 


2.1895 
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a/au 
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图 5 同 图 4, 但 在 引力 模型 中 添加 K = 
0.001 W .m1!.K-! 的 Yarkovsky 效 应 . 


Fig.5 The same as Fig. 4, but add the Yarkovsky effect of 
K = 0.001 W- m~t. KF. 


图 8 是 5 个 在 Yarkovsky 模 型 下 测试 粒子 在 311P 
所 处 共振 区 域 的 停留 时 间 . 纯 引 力 模 型 下 测试 粒 
子 在 50 Myr 内 都 未 脱离 共振 区 域 . 在 Yarkovsky 模 
型 下 , 大 部 分 测试 粒子 在 10 Myr 左 右 脱 离 共振 
区 域 ， 只 有 18 个 测试 粒子 仍 处 于 共振 区 域 . 从 
图 8 第 1 个 bin 中 看 出 , 在 K = 0.001, 0.01 W- m7! 
K-! 模 型 中 测试 粒子 比 其 他 3 个 模型 离开 共振 
区 域 更 快 . 图 中 所 标识 的 中 值 也 可 以 看 出 随 
着 热 导 系数 的 增 大 , 分 布 的 中 值 也 随 之 变 大 ， 
MK = 0.001 W - m~t. K-! 时 的 2 Myr 左 右 到 K = 
0.01 W - m~! .K-71 时 的 3 Myr 左 右 , 再 到 KK = 0.1、 
1.0、5.0 W .m-1.K-! 时 的 5 Myr 左 右 . 


3.2 ”与 大 行星 密 近 交汇 以 及 撞击 的 统计 

本 节 统 计 了 311P 的 500 个 测试 粒子 在 50 Myr 年 
内 与 大 行星 的 撞击 以 及 密 近 交汇 频率 . 

311P 的 500 个 测试 粒子 在 不 同 热 导 系数 的 Yar- 
sells Leg 近 交 汇 统 计 如 网 9. 4k = 
0.001 W- m~t . 开 - 时， 在 20-25 Myr 内 出 现 与 火 
星 密 近 交汇 的 测试 粒子 , 与 火星 密 近 交汇 次 数 随 
时 间 增 加 , 27 个 测试 粒子 a 近 交 汇 1926 次 . 
在 50 Myr 左 右 与 地 月 系统 密 近 交汇 次 数 显 著 增 
大 , 8 个 测试 粒子 与 地 月 系统 密 近 交汇 1160 次 . 通 
过 对 密 近 交汇 次 数 与 Hl 半 径 关系 的 统计 中 可 以 
发 现 , 在 50 Myr 内 311P 进 入 地 月 系统 与 火星 Hil 半 
径 的 次 数 相 当 . K = 5.0 W- m! KSK = 
0.001 W - m-1 :天 -! 情 况 下 的 Hil 半 径 统计 图 的 趋 
势 相 似 , 测试 粒子 和 地 月 系统 的 密 近 交汇 次 数 与 
火星 相当 , EK = 0.001 W . m-1.K-! 时 与 火星 密 
近 交 汇 次 数 略 多 , 而 K = 5.0 W- m. kK WS 
地 月 系统 密 近 交汇 次 数 略 多 . 这 可 能 是 由 于 K = 
5.0 W -m~t .K-! 情 况 下 Yarkovsky 效 应 对 半 长 径 
偏 移 量 的 影响 大 多 为 减 小 . 值得 注意 的 是 , 尽管 与 
火星 和 地 月 系统 总 密 近 交汇 次 数 相 当 , 但 与 地 月 
系统 密 an are 少 于 火星 . K = 
0.001 W .mr- -+ 情况 下 有 8 个 测试 粒子 与 地 月 
系统 密 近 交汇 ， A = 5.0 W.m-l.K-I 情 况 下 仅 
有 一 个 测试 粒子 与 地 球 密 近 交汇 , 且 到 达 共 轨 区 
域 九 :;( 即 与 地 球 1:1 轨 道 共 振 区 域 , 0.994 au < a < 
1.006 au), 因此 311P 变 为 近 地 小 行星 的 概率 虽然 很 
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小 但 不 为 0. 5K = 0.001 W .:m-1.K-! 模 型 相 比 ， 


其 他 4 个 模型 中 测试 粒子 与 火星 密 近 交汇 的 次 数 明 
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16 纯 引力 模型 下 测试 粒子 轨道 经 50 Myr 演 化 在 后 (a, e) ( 左 图 ) 和 (a, i) ( 右 图 ) 空 间 中 的 分 布 情况 , 其 中 (a, e) 平 面 内 , 红色 曲线 为 近 地 天 体 区 域 


边界 gq = 1.3 au, 黑色 实 线 为 g = dearth HQ = Qearth, 黑 色 点 线 为 q = dmars FUQ = Qmars; 黑色 虚线 为 q = Qvenus UQ = Qvenuss 其 中 Q 为 远 


Fig.6 Orbital distribution of the orbital clones in the plane (a, e) 


EB, q., Q. AHN AEM 7 BE H ee Alize H HE. 


left) and (a, i) (right) in the purely gravitational model after 


50 Myr. In the (a, e) plot, the red curve is the boundary of the near-Earth object (NEO) region q = 1.3 au, the black curves are 


q = Qearth and Q = Qearth, the dotted black curves are q = qmars and Q = Qmars and the dashed black curves are q = qvenus and 


Q = Qvenus, where Q is the aphelion distance, and q., Q. are the perihelion and aphelion distances of planets, respectively. 
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) 和 (a, i) ( 右 图 ) 平 面 内 的 轨道 分 布 


417 考虑 K = 0.001 W-m—!-K7! Yarkovsky 效 应 的 引力 模型 下 测试 粒子 经 50 Myr 演 化 后 在 (a, e) (AA 


线 含义 同 图 6. 


Fig.7 Orbital distribution of the orbital clones in the plane (a, e) (left) and (a, i) (right) in the gravitational model considering 


Yarkovsky effect of K = 0.001 W-m~!- K~! after 50 Myr. In the 
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a, e) plot, the meaning of the curves are the same as Fig. 6. 
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48 考虑 Yarkovsky 效 应 的 引力 模型 下 311P 测 试 粒子 共振 区 域 停留 
村 间 分 布 , 相应 的 热 导 系数 在 图 例 中 给 出 . 垂直 虚线 表示 分 布 的 中 位 数 

以 及 对 应 的 数值 在 图 例 中 给 出 . 


Fig.8 Distributions of time spent in the resonance region of 
311P orbital clones under the Yarkovsky effect included in 
the gravitational model, and the corresponding value of 
conductivity is presented in the label. Vertical dashed lines 
indicate the median of the distributions, and the 


corresponding numerical values are also shown in the label. 


表 3 给 出 了 测试 粒子 撞击 、 抛 射 百 分 比 统计 
结果 , 大 多 数 测试 粒子 在 50 Myr 内 存活 下 来 , 但 
有 7% 的 测试 粒子 与 太阳 相 撞 , 其 中 5% 都 是 导热 
系数 KK = 0.001 W . m-!.K-! 情 况 下 的 测试 粒子 . 
2% 的 测试 粒子 由 于 与 木星 的 密 近 交汇 而 从 太阳 系 
抛射 出 去 . 尽管 311P 与 大 行星 的 密 近 交汇 较为 频 
繁 , 但 与 大 行星 碰撞 的 概率 很 小 , 在 500 个 测试 粒 
子 中 没有 与 大 行星 碰撞 而 终止 轨道 演化 的 测试 粒 
子 . 这 是 由 于 对 于 类 地 行星 的 密 近 交汇 不 足以 将 小 
天 体 置 于 双 曲 轨道 中 , 而 该 时 标 内 的 Yarkovsky 效 
应 不 足以 使 得 311P 与 木星 、 土 星 或 海王 星 密 近 交 
汇 而 抛射 至 太阳 系 外 . 表 4 给 出 了 测试 粒子 与 大 
行星 密 近 交汇 的 百分比 , 311P 密 近 交 汇 最 为 频繁 
的 是 火星 , 尤其 是 KK = 0.001 W- m! .KK-I 时 的 
Yarkovsky 模 型 下 有 27% 的 测试 粒子 与 火星 发 生 
密 近 交汇 , 成 为 越 火 小 行星 . 这 一 模拟 结果 与 
Vokrouhlicky 等 8B 的 模拟 结果 一 致 , 对 于 100 m 量 


级 大 小 的 小 天 体 如 果 其 表面 履 盖 一 层 风 化 层 , WE 
要 受 昼夜 效应 的 影响 半 长 轴 可 随机 游 走 至 wm 共 振 ， 
从 而 成 为 越 火 小 行星 或 近 地 小 行星 . 

综 上 所 述 , 311P 将 在 大 约 20-40 Myr 时 间 尺 度 
范围 内 成 为 越 火 小 行星 , 而 成 为 近 地 小 行星 的 概率 
较 小 , 500 个 测试 粒子 中 仅 有 15 个 到 达 过 近 地 区 域 
(q < 1.3 au), 其 中 7 个 测试 粒子 与 太阳 相 撞 , 还 有 
一 个 测试 粒子 到 达 了 地 球 共 轨 区 域 T.1. 


#3 与 太阳 、 火 星 碰 撞 或 抛射 出 太阳 系 的 测试 粒子 百 分 

比 . 最 后 一 列 是 在 50 Myr 积 分 时 标 内 存活 粒子 的 百分比 

Table 3 Percentage of clones that ended their 
evolution by colliding with the Sun, Mars, or by 
being ejected from the Solar System. The last 

column indicates the percentage of clones that 

survived the entire 50 Myr 
Coll. Coll. Other 


K/ Survived 
with with removed 


Oe) Sun Mars objects et 
0.001 5% 0% 1% 94% 
0.01 1% 0% 0% 99% 
0.1 1% 0% 0% 99% 
1 0% 0% 0% 100% 
5 0% 1% 0% 99% 


#4 考虑 不 同 导热 系数 的 Yarkovsky 效 应 下 311P 测 试 粒子 
与 3 个 类 地 行星 金星 、 地 球 和 火星 密 近 交汇 的 频率 
Table 4 Close encounter frequencies between 
311P orbital clones and three terrestrial planets 
(Venus, Earth, and Mars) under Yarkovsky 


effect with varying thermal conductivities 


K/(W.m!.K ') Venus Earth Mars 


0.001 7% 8% 27% 
0.01 2% 2% 16% 
0.1 1% 1% 6% 

1 0% 0% 3% 
5 1% 1% 10% 
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图 9 在 考虑 不 同 导 热 系 数 的 Yarkovsky 效 应 的 引力 模型 下 , 左 列 为 500 个 测试 粒子 与 4 个 大 行星 在 50 Myr 时 标 内 的 密 近 交汇 次 数 统计 图 , 右 列 为 密 


近 交 汇 时 落 在 若干 倍 Hill 半 径 内 的 次 数 统计 图 , 其 中 没有 显示 的 大 行星 均 为 没有 密 近 交汇 . 


Fig.9 Under the gravitational models considering Yarkovsky effect with varying thermal conductivity, the left panels display 
statistics regarding the number of close encounters between 500 orbital clones and four planets within 50 Myr. Meanwhile, the 
right panels depict statistics on the number of encounters occurring within several times the Hill radius of the planets. Planets 


not shown indicate no close encounters. 
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图 9 续 
Fig.9 Continued 
A 讨论 应 的 偏 移 量 为 常量 . 另外 一 个 影响 轨道 偏 移 量 的 
‘ 3 ; ON 为 轨道 偏心 率 . itale 等 6 发 现 对 于 1 量 
本 文 考虑 Yrkvosley 效 应 的 方式 是 在 负 首 速度 eean rig le ee 
‘ g 、 小 的 小 x i ees ARS KT 
方向 施加 了 恒定 的 等 效力 , 这 一 简化 方式 在 轨道 和 级 大 小 的 小 行星 , 当 轨 道 偏 心率 由 0.1 变 为 0.9 时 ， 


自转 轴 演 化 中 有 一 定 的 局 限 性 . nent 企 其 他 因素 
对 半 长 径 偏 移 量 的 影响 . 

在 轨道 演化 方面 , Yarkovsky 效 应 主要 与 日 心 
距 有 关 , 即 受 轨道 的 半 长 径 与 偏心 率 的 影响 29], 当 
天 体 与 太阳 相距 越 近 该 效应 越 强 , 因此 随 着 轨道 
演化 , 半 长 径 偏 移 量 应 该 也 随 之 变化 . 然而 从 上 
述 模拟 中 看 出 Yarkovsky 效 应 主要 影响 主 带 中 位 于 
弱 三 体 共振 处 (a € [2.1874, 2.1912] au) 的 测试 粒子 ， 
当 其 成 为 越 火 小 行星 或 近 地 小 行星 时 , 引力 效应 
为 主要 影响 半 长 径 的 因素 , 因此 假设 Yarkovsky 效 


XR 
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Yarkovsky 昼 夜 效 应 会 增 大 一 个 量 级 左右 . 此 时 恒 
定 偏 移 量 可 能 不 再 适用 , 但 图 7 说 明 仅 有 一 小 部 分 
( 约 1.4%) 与 太阳 相 撞 和 抛射 出 太阳 系 的 测试 粒子 
(如 表 5) 达 到 高 偏心 率 并 结束 轨道 演化 . 同时 图 7 也 
说 明了 测试 粒子 在 共振 带 中 偏心 率 大 都 小 于 0.3， 
因此 由 于 高 偏心 率 导 致 Yarkovsky 偏 移 量 数量 级 的 
变化 并 不 会 导致 结果 的 改变 , 尤其 是 对 于 测试 粒子 
在 共振 带 区 域 的 演化 . 

在 自转 轴 的 演化 方面 YORP 效 应 Rl 291, 4E 
破坏 性 碰撞 B93 以 及 密 近 交汇 9 等 都 会 改变 小 行 
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星 的 旋转 状态 . 未 遭受 碰撞 的 小 行星 一 般 主要 受 
YORP 效 应 的 影响 , 一 旦 经 历 非 破坏 性 碰撞 则 会 显 
著 改 变 自转 参数 . StatlerB49 关 于 YORP 效 应 对 小 尺 

度 天 体 表 面 特征 分 析 表 明 自 转轴 演化 很 可 能 是 


能 是 随 
机 的 , 并 且 之 后 Cotto-Figueroa 等 B5 也 通过 硬 球 离 
散 元 数值 


模拟 (hard-sphere discret element numer- 
ical simulation) 证 实 了 这 一 观点 . WiegertB3 同 样 
认为 非 破坏 性 磁 撞 会 重新 随机 定向 自转 轴 , ALE 
非 破 坏 性 碰撞 与 YORP 效 应 同样 重要 , 尤其 是 对 于 
像 311P 这 种 小 尺度 的 小 行星 . 根据 Farinella 等 B9 关 
于 撞击 频率 的 研究 , 小 行星 在 主 带 中 由 于 碰撞 导致 
自转 轴 随 机 定向 的 时 标 与 小 行星 的 大 小 和 自转 周 
期 有 关 , 311P 这 种 尺寸 大 小 的 天 体 在 100 Myr 的 时 
标 下 不 会 频繁 发 生 撞 击 事 件 , 即 不 会 频繁 随机 定向 
自转 轴 方 向 , 这 远大 于 模拟 中 测试 粒子 脱离 共振 区 
域 的 时 标 . 因此 辜 撞 以 及 密 近 交汇 在 演化 中 影响 自 
转轴 分 布 的 程度 有 限 . 

而 YORP 效 应 的 估计 则 要 复杂 得 多 , 因此 对 
其 做 了 进一步 的 模拟 与 讨论 . YORP 效 应 是 将 倾 
角 7 推 向 渐 近 状态 之 一 (0"、90? 或 180?) 的 热 扭矩 ， 
同时 会 导致 小 行星 旋转 速率 减缓 或 加 快 . 然而 这 一 
效应 与 天 体 的 形状 和 热力 学 性 质 密切 相关 , 对 于 
大 部 分 小 天 体 这 些 参数 都 是 未 知 的 . 因此 本 文采 
Fenucci 等 85 的 静态 YORP 模 型 , 小 天 体 倾角 y 和 
自转 速率 w' = 1/P (由 于 近日 点 幅 角 使 用 ww 表示 ， 
这 里 用 w' 表 示 自 转速 率 ) 的 演化 由 下 式 给 出 : 


U 


— 


I =f Js 
T 7 (5) 
dy _ 9(y) 
ta~ w’ 


其 中 f、g 为 扭矩 大 小 , f 和 g 的 选择 可 参考 Fenucci 
4015] YORP 效 应 通过 影响 Yarkovsky 效 应 来 影响 
轨道 动力 学 演化 因此 轨道 演化 过 程 中 , 在 每 个 积 

分 步 长 处 都 会 通过 (5) 式 更 新 y 和 PP 来 更 新 (2) 式 中 
的 半 长 径 偏 移 量 . 轨道 演化 的 初始 条 件 设置 与 2.2 
节 所 述 一 致 . 由 于 KK = 0.001 W .m71.K-! 为 导 
臻 半 长 径 偏 移 量 最 大 的 情况 , 这 里 只 评估 该 情况 
下 的 Yarkovsky 模 型 与 YORP 模 型 的 差别 . 311P 在 
Yarkovsky 模 型 与 YORP 模 型 下 的 共振 区 内 时 标 
分 布 如 图 10, 该 分 布 可 很 好 地 拟 合 为 对 数 正 态 分 


N 


m| 
ZN 
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ti, BUX = et, 其 中 2 为 标准 正 态 分 布 . Yar- 
kovsky 模 型 下 的 分 布 参数 为 4/ = 一 0.689、o = 
0.830. YORP 模 型 下 的 分 布 参数 为 / = 一 0.650、 


o = 0.806. 因此 两 个 模型 下 测试 粒子 离开 311P 共 
振 区 域 的 时 标 分 布 几乎 一 致 . 


2 
w 
© 


Probability Density Function 


Probability Density Function 


20 


10 


15 
t/Myr 
到 分 别 为 六 = 0.001 W-m!. KIN 
# 振 区 域 停留 时 


10 上 图 和 下 
虑 Yarkovsky 效 应 和 YORP 效 应 的 引力 模型 下 311P 
标的 对 数 正 态 分 布 拟 合 情 况 . 


Fig.10 The top and bottom panels show the fit of the 


log-normal distribution of resident time in the 311P 
resonance region under the gravitational models considering 
the Yarkovsky effect and YORP effect for the case 
K = 0.001 W- m7! - KT! respectively. 


另 一 方面 Golubov 等 B14 指 出 小 天 体 在 经 历 


几 次 YORP 循 环 后 , 其 自 旋 状态 会 达到 一 个 平衡 
点 , 即 渐 近 态 . 为 了 评估 311P 自 转轴 到 达 渐 近 态 


以 及 自转 周期 变 为 自转 周期 分 裂 极限 的 时 标 ,对 
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测试 粒子 进一步 做 了 10 Myr 的 轨道 演化 . 其 他 初 
始 条 件 保 持 不 变 , 初始 自转 周期 设 为 已 = 6h, 轨 
道 演 化 时 间 10 Myr. 结果 如 表 5 与 图 11 所 示 . 自转 
周期 由 6 h 演 化 至 1.6 h 的 时 标 中 值 约 为 2 Myr. 由 
于 对 311P 的 模拟 显示 到 达 渐 近 态 的 时 间 大 致 为 演 
化 的 前 1 Myr 左 右 , 半 长 径 的 偏 移 量 大 小 不 超过 
0.02 au, 因此 YORP 效 应 对 半 长 径 的 整体 偏 移 
响 较 小 , 不 会 对 模拟 的 最 终结 果 造 成 太 大 影响 . 


=E E 


E 


表 5 考虑 不 同 导热 系数 的 YORP 效 应 的 引力 模型 下 测试 粒 


子 达到 渐 近 态 与 自转 周期 分 裂 极限 的 概率 与 时 标 中 值 统计 
Table 5 The probability and median timescale 


statistics of orbital clones reaching the 
asymptotic state and spin-period splitting limit 
under the gravitational model considering the 
YORP effect with different thermal conductivity 


K/(W-m KR) tga Piens trot? 
0.001 950 kyr 94% 2.03 Myr 
0.01 1Myr 88% 1.83 Myr 
0.1 691 kyr 82% 2.65 Myr 
1 906 kyr 89% 1.80 Myr 
5 1.3 Myr 94% 2.41 Myr 


* Thermal conductivity; 

> The median timescale for orbital clones to reach 
the asymptotic state; 

© The probability of orbital clones reaching the spin- 
period splitting limit of 1.6 hours; 

4 The median timescale for orbital clones to reach 


the spin-period splitting limit. 


结论 
本 文通 过 对 500 个 测试 粒子 向 前 积分 50 Myr 研 
究 了 311P 的 长 期 动力 学 演化 , 并 且 考 虑 了 Yarkov- 
sky 效 应 对 其 半 长 径 演 化 的 影响 . 模拟 结果 如 下 : 

1. Yarkovsky 效 应 仅 在 低 导 热 系 数 下 对 其 动 
力学 稳定 性 造成 明显 影响 , 与 纯 引 力 模型 相 比 ， 
Yarkovsky 效 应 加 速 了 311P 离 开 共振 带 , 在 5 Myr 
时 标 下 约 有 一 半 测 试 粒子 离开 共振 带 , 在 10 Myr 时 
标 下 有 约 80% 左 右 测试 粒子 离开 共振 带 ; 


5 
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K=0.001 W- m-1. K7} 
K=0.01W-m?-K? | 
K=0.1 W- m-1. K-! 
K=1.0 W- m~}. K-! 
K=5.0 W: m~}. K-! 
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Number of objects 
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311P 测 试 粒子 在 不 同 热 导 系数 的 YORP 模 型 下 达到 
分 裂 极 限时 标的 统计 


11 转 周期 


R] 


Fig.11 The timescale statistics of 311P orbital clones 
reaching the spin-period splitting limit under the YORP 


model with different thermal conductivity 


2. 有 12.4% 的 概率 在 20-40 Myr 时 间 内 通过 zx 
共振 , 可 能 成 为 越 火 小 行星 , 但 在 50 Myr 内 成 为 近 
地 小 行星 的 概率 较 小 , 到 达 地 球 共 轨 区 域 的 概率 几 
FXO; 

3. 在 低热 导 率 (覆盖 风化 层 的 岩 质 小 行星 ) 的 
情况 下 , 311P 受 Yarkovsky 昼 夜 效应 的 影响 成 为 越 
火 小 行星 的 概率 更 大 , 约 为 27%; 

4. 如 果 YORP 效 应 是 311P 旋 转 不 稳定 的 原因 ， 
其 旋转 裂变 的 时 标 大 约 为 2 Myr; 


=) 


那么 


5. 密 近 交汇 、 非 破坏 性 碰撞 和 YORP 效 应 对 
311P 轨 道 演 化 的 影响 较 小 , 因此 在 50 Myr 时 标 内 


可 只 考虑 Yarkovsky 效 应 对 其 动力 学 的 影响 . 

综 上 , 在 各 种 因素 的 综合 考虑 下 ，311P 将 在 
10 Myr 时 标 内 离开 共振 带 , 同时 在 表面 覆盖 风化 层 
的 情况 下 可 能 成 为 越 火 小 行星 , 成 为 近 地 小 行星 的 
概率 不 大 . 
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Dynamical Characteristics of Active Asteroid 


311P/PANSTARRS 
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ABSTRACT Active asteroids have received lots of attention since they were discovered. 311P/PANSTARRS 


is one of the most interesting active asteroids in China recently and one of the targets of Tianwen-2. Because 
of its small size of about 400 m, the Yarkovsky effect may have a significantly influence on its long-term 
dynamics. This paper discussed the changes in the long-term motion of 311P/PANSTARRS caused by the 
Yarkovsky effect. By assuming different surface compositions, this simulation introduced the semi-major 
axis drift by propagating orbits of orbital clones, and the effects of other factors such as close encounters, 
meteoroid impacts and YORP effect on the orbit evolution of 311P/PANSTARRS were also discussed. 
Additionally, we estimated the time scale for 311P/PANSTARRS to reach its rotation period of splitting 
limit. The results of the simulations show that the Yarkovsky effect may cause 311P/PANSTARRS to 
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exit from the resonance region faster when compared to a purely gravitational model. 311P/PANSTARRS 
will leave the current resonance region after roughly 10 Myr and have a chance to become a Mars-crossing 
asteroid through ve secular resonance due to the diurnal Yarkovsky effect if its surface is covered by a re- 
golith layer. It is concluded that 311P/PANSTARRS is stable at least 10 Myr time scale even if taking the 
Yarkovsky effect and the YORP effect into account. Furthermore, the YORP effect may not significantly 


affect the semi-major axis drift of 311P/PANSTARRS. 


Key words celestial mechanics, methods: statistical, planets and satellites: dynamical evolution and 
stability, minor planets, asteroids: individual: 311P/PANSTARRS 


附录 


本 文 利 用 Yarkovsky 效 应 的 昼夜 效应 和 季节 效 

应 来 计算 相应 的 半 长 径 变 化 率 , 并 对 其 进行 蒙特 卡 
罗 实 验 . 昼夜 效应 和 季节 效应 导致 的 半 长 径 变化 率 
由 下 式 给 出 [3]; 

da 8a P j 

( =- F(R 0) cosy, (6) 

da\ 4a © , nt 

ee = va Coo y, (7) 


表面 的 吸收 系数 , OAK RM, © 
RW. JW 表示 距离 太阳 a (au) 处 小 行星 所 收 到 
的 太阳 辐射 流量 ,cv 为 真空 中 的 光速 ,Rs 和 ma 分别 
为 小 行星 的 半径 和 质量 . ce, 为 平 运动 角速度 ，? 为 
自转 轴 倾 角 , 下 是 关于 小 行星 半径 和 热 物理 参数 相 
关 的 函数 . 丸和 忆 分 别 为 小 行星 半径 缩放 的 无 
HEL, 定义 为: 


其 中 , a 为 


_ xR2W, 


量 
Æ, Cv 


里 


R R 
Ra = 天 ,及 = 于， (8) 
la ls 
FEAL Ag DA AIA PL IM A EY IN PBS BY PIE 
深度 , BẸ: 
K K 
[= 9 
a pC Wrot : PCWrev ( ) 
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1 和 1 由 天 体 表 面 热 导 率 玉 、 比 热 容 C 和 小 行星 密 
度 p 确 定 . 其 中 自转 频率 wo 用 于 昼夜 效应 , wev 用 
于 季节 效应 . 小 行星 热 物 理 参 数 相 关 量 6 和 @。 由 
下 式 决 定 : 
VPK Cwrot VPK Cwrev 
Ga 一 , Os = . (10) 
Eopls3 EOb73 


ob 为 Stefan-Boltzmann 常 数 , = 为 表面 发 射 率 , TN 
日 下 点 温度 , HeoT* = aWi 计 算 . (6) 式 和 (7) 式 中 
的 取决 于 相应 的 半径 缩放 量 和 热 物 理 参 数 , 并 且 
决定 了 半 长 径 偏 移 的 量 级 . 的 表达 式 为 : 


F 


kı (R) © 


EOFS 1+2 k (R) O + ks 


mya O” 


Rikon ka kst N Aa IE 42 H E fT R A, 它 
们 的 定义 在 Vokrouhlickyl8s-39 中 已 详细 给 出 . 注 
意 ， AFPRNAK, 季节 效应 总 是 使 得 半 长 径 减 
小 , 且 在 倾角 达到 90°* 时 达到 最 大 , 在 0° 或 180°? 时 增 
加 量 为 0. 而 昼夜 效应 则 在 y < 90? 时 半 长 径 增 大 ， 
在 7 > 90° 时 半 长 径 减 小 , 在 倾角 为 0? 或 180? 该 效应 
作用 最 大 , 在 倾角 为 90" 时 增加 量 为 0. 
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